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Abstract 

For a dissipative system with Ohmic friction, we obtain a simple and exact solution 
for the wave function of the system plus the bath. It is described by the direct 
product in two independent Hilbert space. One of them is described by an effective 
Hamiltonian, the other represents the effect of the bath, i.e., the Brownian motion, 
thus clarifying the structure of the wave function of the system whose energy is 
dissipated by its interaction with the bath. No path integral technology is needed 
in this treatment. The derivation of the Weisskopf-Wigner line width theory follows 
easily. 
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The simplest example of a dissipative system, an harmonic oscillator coupled to 
the environment, which is a bath of harmonic oscillators, has been the subject of 
extensive studies^"^^. We shall show in the present paper that in a special case, the 
Ohmic case (to be defined later), the dissipative system can be exactly treated both 
classically and quantum mechanically, thereby clarifying the sense in which the wave 
function is describable by an effective Hamiltonian. In this treatment path integral 
technology is not needed, and our presentation is self-contained. 

We consider the problem discussed by Caldeira and Leggett (CL)^, an harmonic 
oscillator system (the dissipative system) with coordinate q, mass M, and frequency 
{uq + Aa;^)^/^, interacting with a bath of N harmonic oscillators of coordinates Xj, 
mass rrij, and frequency utj., where Aa;^ is a shift induced by the coupling already 
discussed by CL. The Hamiltonian of the system and the bath is: 

The dynamic equation for operators in the Heisenberg representation leads to the 
following set of equations of motions: 

Mq = -MuJlq + MAu?q - ^ CjX^ , (2) 

3 

rrijX — —mjUj'^Xj — cjq {j — 1, 2, ....N). (3) 

Now, applying the Laplace transform^ (the bars are our notations for the Laplace 
transform, s is the Laplace transform of time t), equations (2), (3) can be used to 
eliminate bath variables xj to obtain the equation for q 

M{s-q-sq,-q,)^-Mulq-MAu-q-Y.c,-^^^f^-^ ^ J^^, (4) 

3 3 j J \ J / 

where qo, qo, Xjo, Xjo are the initial values of the respective operators in the Heisen- 
berg representation. Assuming the number of bath oscillators is large enough so that 
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we can replace the sum over j by an integration over ujj, the coefficient of the last 
term can then be separated into two terms: 



Jo rUjUJj Jo mjUJj{s'^ + ujj) ' 



where p{u!j) is the bath oscillator density. Following an argument similar to the one 
pointed out by CL^, the requirement that the system becomes a damped oscillator 
with frequency cuq and damping rate r) in the classical limit, known as the "Ohmic 
friction" condition, leads to the following constraint: 



2r]Mmjuj 

Pi^j) = — 3^- (6) 



TT 

By observing eq.(4) and eq.(5), it can be shown that the second term of eq.(5) leads 
to damping with the damping constant r], while the first term of eq. (5) represents a 
frequency shift. If the frequency renormalization constant Ao;^ is chosen to satisfy: 

MA.= =r^M!i, (7) 

the frequency is shifted to luq. Then equation (4) is simplified, and its inverse Laplace 
transform gives the quantum Langevin equation, valid at time t >= 0+: 

q + m + ulq^ f{t) , (8) 

with the Brownian motion driving force: 

f{t) = -Y, Cj {xjo cosujj + XjQ ^ ) . (9) 

During the derivation, in order to carry out the integral in eq.(5), we used the re- 
quirement of the inverse Laplace transform that s must pass all the singular points 
from right of the complex plane, and hence real(s) >0. 

Equation (8) and (9) are the equations of a driven damped harmonic oscillator, 
the solution of which is well known as a linear combination of the initial values at go, 



qo, Xjo, and Xjo : 

q{t) = ai{t)qo + 02(^)^0 + J2 (^ji(*)^jO + bj2{t)xjo) , (10) 

3 

Xi{t) = aii{t)qo + ai2{t)qo + '^{Piji{t)xjo + (3ij2{t)xjo) , with (11) 

i 

ai W = , «2W = , 12 

V — \x u — II 

H = - + uji , and v — - — cu i , (13) 

and here cu — {cUq — rf /AY/"^ is the frequency shifted by damping. (All formulae are 
correct whether uj is real or imaginary. To avoid a minor detail of the initial value 
problem, we have redefined the initial time as t=0+.) The explicit expression for hji, 
bj2, CKii, ai2, Piji, and Pij2 are well known in freshmen physics. 

We emphasize that the use of the Laplace transform instead of the Fourier trans- 
form allows us to express q(t) and Xj{t) explicitly in terms of the initial values, as in 
eq.(lO) and eq.(ll). 

The equations (10) and (11) serve as the starting point of subsequent discussions. 
We will proceed to find the Green's function of the full system, and hence the solu- 
tion of the wave function in Schoedinger representation. The result tells us in what 
sense the damped oscillator is described by an effective Hamiltonian without the bath 
variables and gives its specific form, it also shows that under this condition, the wave 
function can be factorized, and the main factor relevant to the damped oscillator is 
a solution of the Schoedinger equation with an effective Hamiltonian. 

Equations (10) and (11) are correct both in classical mechanics and in quantum 
mechanics in the Heisenberg representation. We notice that q(t) and Xj{t) are both 
linear superpositions of qo, %, Xjq, Xjo with c-number coefficients. The commutation 
rules between q(t), q{t), Xj{t), Xj(t) are [q{t),q{t)] = [xj{t), Xj{t)] = ^, and 
operators q(t) and q{t) commute with Xj{t) and Xj{t). One can prove these commu- 
tation rules by two ways: (a). By direct computation, using the fact that at t=0, 
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they are correct, (b). By the general principle that q{t), q{t), Xj{t), Xj{t) are related 
by a unitary transformation to Qq, Qq, Xjo, andxjo. 

Equation (10) and (11) show that the operators q(t) and Xj{t) can each be written 
as a sum of two terms: 

= W + EOW , (14) 

j 

Xi{t) = Q{t) + Y.Xij{t) , (15) 
j 

where Q{t) and Q{t) are linear in go and go and independent of Xjq and Xjq, and and 
Xij{t) are linear in xjo and Xjq, and independent of go and go. Thus Q(t) and (i{t) 
are operators in one Hilbert space Sq, while ^j{t) and Xij{t) are in an independent 
Hilbert space Sx, and the full Hilbert space is a direct product Sq(S)Sx- 
We shall first analyze the structure of Sq. We write that: 

tfi 

Q{t) = aiQo + a2Qo = aiQo - a2 yj^TT ■ W 

To explicitly show that we are discussing the Sq space, we define Qo = 9o, Qo = ?o- 
The eigenfunction of Q(t) with an eigenvalue denoted by Qi, in the Qo representation, 
is easily calculated to be: 



(17) 



2ha2 

with (f) as an arbitrary phase, i.e., a real number. This eigenfunction is related to the 
Green's function G(Qi,Qo;^,0) =< Qi\UQ{t)\Qo >,hy UQ^{Qo,t) =< Qo\UQ\t)\Qi > 
=< Qo\UQ(t)\Qi >= G*{Qi,Qo;t,0), where we denote the evolution operator by 
U{t), which is unitary when we choose the eigenvectors of Q(t) to be orthonormal. 
Thus we have: 



^ iM 



2ha2 



{aiQl-2QoQi + (l)(Qi,t)) 



(18) 
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Next, we shall determine the arbitrary phase ^(Qi,^), which is the phase of the 
eigenvectors of Q(t). Using eq.(16), we find the commutation rule for Q and Q: 

[Q, Q] = (aia2 - aia2)[go, Qo] = e-'^'iQo, Q^] = ^-^ih . (19) 

Thus we define the canonical momentum P(t) as: 

P{t) = Me"'Q{t) = M e''\aiQo + a^Qo) = aiQo - «2^^ , (20) 

and get the commutation rule: [Q{t),P{t)] = iU. The eigenfunction of P(t) can be 
calculated in two ways: (a). We can calculate the eigenvector of P(t) in the Qq 
representation using eq. (20) and then use the Green's function eq.(18) to transform 
it into the Q(t) representation; (b). The commutation rule [Q(t), Pit)] — ih requires 
that P(t)=— i/i^, so the eigenfunction of P(t) with eigenvalue Pi is exp[i^Q]. By 
comparing these two solutions, the arbitrary phase (f){Qi,t) in the Green's function is 
determined to within a phase (f){t), which is independent of Qi. (f){t) is an arbitrary 
real function of time, except that 0(0) = so that it satisfies the condition that at 
t=0, the Green's function becomes S{Qi — Qo)- Thus we obtain the Green's function 
in the Sq space: 



2ha2 ^ ^^^^ h 



(21) 

It is then straight forward to derive the Hamiltonian Hq using the following rela- 
tion: 

Hq = in^-^UQ^ , (22) 

and remembering that the matrix elements of Uq and Uq^ are the Green's function 
and its conjugate. The result is: 



Since is arbitrary except that 0(0) = 0, we can take 0(t) = 0. Therefore we 

have derived the well known effective Hamiltonian for the dissipative system. We 

emphasize that the expression for Hq is here derived, while in usual literature it is 

introduced by more or less heuristic arguments. 

Next, we shall analyze the effect of the bath. Similar to eq. (17) we obtain 

the eigenfunctions 9^.^{xjo,t) for ^j. Using Dirac's notation we have: Q\uq^ >— 

Qi\uq^ >, >= >• Thus |mq > 1% > is eigenvector of q(t), 

j 

with eigenvalue of Q + ^Cj- other words, the eigenvector of q(t) with eigenvalue 

i 

q is \Q, {Cj} >= _ > n 1% >■ 

j 

If initially the wave function is |^o >= iV'o > U \Xjo >, to calculate the wave 

j 

function at time t, we should expand ^'o in terms of the eigenvectors \q, {^j} >, i.e., 
we should calculate 

V^(Q, t) =< uq\iIjo >= j u*Q{qo)ilJo{qo)dqo , (24) 

Xj (Ci, t) =< % Ixio >= j^l {xjo)Xjo{xjo)dxjo ■ (25) 
Then the wave function in the Schoedinger representation at time t is 

^(q,{^j},t)^<q,{^j}\^o>^HQ-EM UxAM ■ (26) 

j j 

Notice that ip{Q, t) of eq.(24) is the wave function in the Schoedinger representa- 
tion with the effective Hamiltonian eq.(23). Hence we have connected the effective 
Hamiltonian approach to the dissipative system problem with the other approaches 
that take both the system and the bath into account. We also notice that even though 
our ^(g, {Cj},t) is in a different representation from that of ^(g, {xj},t), the usual 

probility interpretation is still valid: j j ... J \^{q, {^j}, ^) |^ d^j is the probability 

j 

density of finding the particle at q. Since this solution is very simple, it provides 
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a simple way to analyze other complicated problems, e.g., study the influence of 

Brownian motion on interference, which we shall not elaborate because of space. 

Under certain conditions, eg., at low temperature and when the system q is in 

highly excited states so the range of q is large enough that for all the values of 

which do not have vanishingly small probability, q ^ I X/^jI' ^® approximately 

j 

write: ^{q,{^j},t) = ij{q:t)Y]_Xji^j:'t)- That is, the wave function is factorized, 

j 

the dissipative system q can be described by the wave function ip{q.,t) only, and the 
Brownian motion can be ignored. Therefore it is interesting to examine the width 
of the argument of the wave function V', due to the Brownian motion, i.e., the mean 
value of (X (^jY at time t. It can be calculated using its expression in the Heisenberg 

3 

representation as introduced by (10) and (18). At temperature T, this width is : 

< (zmf >= E ^(\bnit)\' + -]\bAm cot (^) . (27) 

This width is zero initially, then approaches its final equilibrium value in a time 
interval of the order of I/77. At low temperature limit, the equilibrium width is 
simplified to: 

- + arctan ( —] . (28) 
2 V^^/J 

If the damping rate rj is much smaller than the frequency of the oscillator cuq, this 
width happens to become the same as the width of the ground state of the system 
h/ (2ma;). 

Finally, it is interesting to see the distribution of the dissipated energy of the 
harmonic oscillator in the bath, and check if it agrees with the Weisskopf-Wigner line 
width theory^^. For simplicity, we assume zero temperature. To calculate the energy 
dissipated by the system into the j'th bath oscillator, we use eq.(ll) and its derivative 
to obtain the expression of xj (t) and pj = rrijXj (t) , which are then substituted into the 
expression of the energy of the j'th oscillator: hj = + ^rrijLUjX'^. We then calculate 
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a^(t = 00) = 

zTrmw 



the expectation value of hj, assuming initially the system is in the n'th excited state, 
and the bath oscillators in the ground state. We calculate the contribution to the 
expectation value of hj from the system, by keeping only terms which depend on qo, 
and qo- The result is then the energy disspation by the system into the j'th oscillator. 
When multiplied by the density of states p{ujj) eq.(6), it gives the dissipated energy 
spectrum. It is a function of cuj, with a narrow peak near the resonance cuj — uj, if 
the damping is small, i.e., if ?] <^ o;. It is oscillatory with a frequency of 2uj. Its time 
average over a period is : 

Ej = (n + hfi^ 1^-7^ (29) 

2 TT 2_ + _ ^^.)2 

where A^uoj) varies slowly near the resonance, i.e., over the width r]/2 of the peak it 
changes very little. Therefore, we can replace cuj by cu near the peak, and the result 
is simplified to: A{u;) fa 1. Thus eq.(29) shows that the dissipated energy has a 
Lorenzian distribution near the resonance, in agreement with the Weisskopf-Wigner 
line width theory. 
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